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Abstract 

We study the spatial decay of spin polarized hot carrier current in a spin-valve 
structure consisting of a semiconductor quantum wire flanked by half-metallic 
ferromagnetic contacts. The current decays because of D'yakonov-PereF spin 
relaxation in the semiconductor caused by Rashba spin orbit interaction. The 
associated relaxation length is found to decrease with increasing lattice tem- 
perature (in the range 30-77 K) and exhibit a non-monotonic dependence on 
the electric field driving the current. The relaxation lengths are several tens 
of microns which are at least an order of magnitude larger than what has been 
theoretically calculated for two-dimensional structures at comparable temper- 
atures, Rashba interaction strengths and electric fields. This improvement is a 
consequence of one-dimensional carrier confinement that does not necessarily 
suppress carrier scattering, but nevertheless suppresses D'yakonov-Perel' spin 
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Ever since the discovery of the spin- valve effect [1], there has been considerable interest 
in studying spin transport in non-magnetic materials in which spin polarized carriers are 
injected from a ferromagnetic contact and detected by another ferromagnetic contact. The 
spin valve structure has also been employed to devise novel spintronic devices, such as the 
so-called spin field effect transistor [2], where an electron's spin (rather than its charge) is 
employed to elicit transistor action. 

The basic spin valve geometry is shown in the top panel of Fig. f. It consists of a 
semiconductor channel (assumed to be quasi one-dimensional for this study) flanked by 
two half-metallic ferromagnetic contacts. One contact (called the "source") injects spin 
polarized current into the channel and thus acts as a "spin polarizer". The other contact 
acts as a "spin analyzer" and is termed the "drain". Carriers drift from the source to 
the drain under the influence of a driving electric field. When they arrive at the drain, 
they are transmitted with a probability \T\ 2 = cos 2 (9/ 2) where 9 is the angle between 
the electron's spin polarization at the drain end and the drain's magnetization [2]. With 
increasing degree of spin depolarization in the channel (caused by spin relaxation), the 
average "misalignment angle" 9 (for the electron ensemble) increases and consequently the 
transmitted current decreases. Ultimately, when there is no residual spin polarization in the 
current (i.e. carriers are equally likely to have their spins aligned parallel or anti-parallel to 
the drain's magnetization), the transmitted current will fall to 50% of its maximum value. We 
are interested in finding how the (transmitted) spin polarized current falls off with distance 
along the channel at different driving electric fields and temperatures. 

Spins depolarize in the channel primarily because of spin-orbit interactions caused by bulk 
inversion asymmetry (Dresselhaus spin-orbit coupling) [3] and structural inversion asymme- 
try (Rashba spin-orbit coupling) [4]. These spin-orbit couplings are momentum dependent, 
and because different electrons have different momenta that change randomly due to scat- 
tering, the spins become randomized by scattering and the ensemble averaged spin and spin 
polarized current decay with distance. This mechanism of spin relaxation is the D'yakonov- 
Perel' mechanism [5] which is overwhelmingly dominant in quasi one- dimensional structures 
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over the Elliott- Yafet [6] or Bir-Aronov-Pikus [7] mechanisms. The spatial decay of spin due 
to D'yakonov-Perel' mechanism was studied in the past by Bournel, et. al. [11] and Saikin, 
et. al. [12] in two-dimensional channels. They mostly dealt with low driving electric fields 
so that transport is linear or quasi-linear. In contrast, we have studied the spatial decay 
in quasi one-dimensional structures of both spin and spin polarized current at high driving 
electric fields of 1-10 kV/cm, which result in hot carrier transport and non-linear effects. 

In a one-dimensional structure, the spin polarized current due to one electron is pro- 
portional to qv x \T\ 2 where v x is the ensemble averaged velocity of the electrons along the 
channel. As stated before, the quantity \T\ 2 depends on the component of the electron's 
spin polarization along the magnetization of the drain. We will assume that the source and 
drain are both magnetized along the channel's axis (x-axis). This results in the initial spin 
orientation to be along the channel axis. Accordingly, 

|T| 2 = cos 2 (9/2) 

cos(9) = S x /y/S* + S* + S* = S x , (1) 

where S n is the spin component along the n-axis and S x is the normalized value of S x . 
The ensemble averaged spin polarized current at any position x is given by 

I s (x) =qJ2_ f(v x ,S x ,x)v x \T(S x )\ 2 , (2) 

where the velocity (v x )- and spin (S x )- dependent distribution function f(v x ,S x ,x) at any 
position x is found directly from the Monte Carlo simulator described in ref. [13] (all perti- 
nent details of the simulator can be found in ref. [13] and will not be repeated here). We only 
mention that in the simulator, we use a parabolic energy versus velocity dispersion relation 
E = (h 2 /2m*)(nn/W z ) 2 + (h 2 /2m*)(n /W y f + (l/2)m*v 2 (n is the subband index in the 
z-direction) , neglecting any band structure non-parabolicity which is not important in the 
energy ranges encountered [8]. This dispersion relation allows us to calculate the velocity 
v x from the carrier energy E and subband index n (which are tracked in the simulator) 
very easily. If instead we used the energy versus wavevector relation (which is traditional) 
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and then attempted to find v x from the velocity versus wavevector relation, it would have 
been immensely complicated. The reason is that the velocity (or energy) versus wavevector 
relation is spin- dependent in the presence of the Rashba effect [9] and becomes even more 
complicated if the Rashba effect is strong which leads to spin mixing effects [10]. These 
complications would have been overwhelming in our case since we have a continuous dis- 
tribution of spin and hence would have been faced with a denumerably infinite number of 
energy versus wavevector relations. The way to avoid this daunting complication (and the 
associated numerical cost) is to use the energy-velocity relation which is spin-independent 
instead of the energy-wavevector relation which is spin- dependent. 

In the simulation, carriers are injected into a quasi one-dimensional GaAs channel of 
rectangular cross section 30 nm x 4 nm (see top panel of Fig. 1). We have assumed 
that there is a transverse electric field of 100 kV/cm (in the y-direction) that gives rise to 
a structural inversion asymmetry and induces a Rashba effect in the channel. This field 
perturbs the subband energies in the channel but only slightly. The transverse voltage drop 
over a 4 nm wide channel due to this field is 40 meV while the lowest subband energy is 355 
meV. Therefore, the perturbation is 11% for the lowest subband and progressively decreases 
for higher subbands. Consequently, we neglect this perturbation. Electrons are injected 
from a Fermi Dirac distribution with their spins all aligned along the channel axis (x-axis) 
in order to simulate the spin polarizer. At any given position x, we find the spin vector 
S x and compute the quantity \T(S X )\ 2 for every electron using Equation (1). We also find 
the velocity v x for every electron at position x and then compute the spin polarized current 
I s by performing the ensemble averaging given by Equation (2). We have found I s versus 
position x for four different channel electric fields of 1, 2, 4 and 10 kV/cm and two different 
temperatures of 30 and 77 K. 

In the Monte Carlo simulation, we find that typically the lowest 4 subbands are occupied 
by electrons. The Dresselhaus interaction strength is different in different subbands because 
the interaction strength is proportional to (nir/W) 2 where W is the transverse dimension 
(4 nm in our case) and n is the subband index. As electrons are scattered between various 
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subbands, they witness varying spin orbit interaction which results in a decaying envelope of 
the ensemble averaged spin (or spin polarized current). This is the cause of D'yakonov-Perel 
relaxation in quasi 1-d structures. 

In Fig. 1, we show the spatial decay of the normalized spin polarized current I s for the 
four different (x-directed) channel electric fields at a temperature of 30 K. In Fig. 2, we show 
the same quantity (along with the spatial decay of the ensemble averaged spin component 
S x ) at an electric field of 2 kV/cm at temperatures of 30 and 77 K. Spin depolarization is 
complete when I s reaches a value of 0.5. At this point, an electron is equally likely to have its 
spin aligned parallel or anti-parallel to the drain's magnetization (and therefore it is equally 
likely to be transmitted or blocked). We can define a "relaxation length" as the distance over 
which the injected spin polarized current decays to 50% of its value (i.e. becomes completely 
depolarized). Table I gives the relaxation lengths at different electric field strengths and 
different temperatures. 

As expected, the relaxation length decreases with increasing carrier temperature because 
of increased scattering that causes increased spin depolarization. The electric field, on the 
other hand, has two opposing effects. The scattering rate increases slowly with the electric 
field, but so does the ensemble averaged carrier drift velocity until the saturation velocity is 
reached. A larger drift velocity makes the carriers travel a greater distance before getting 
depolarized. Consequently, the relaxation length at first increases with increasing electric 
field, but once the drift velocity begins to saturate, the increased scattering takes over and 
the relaxation length starts to decrease with increasing electric field. The dependence of 
relaxation length on the electric field is therefore non-monotonic. 

Based on the data in Table I, we find that the relaxation length for spin polarized current 
is very large (between 20 and 100 /iin for the cases considered). This is at least an order of 
magnitude larger than what was calculated for two-dimensional structures [11,12] at compa- 
rable temperatures and driving electric fields. This difference is not due to any suppression 
of scattering. In fact, even though elastic scattering is suppressed in quasi one-dimensional 
structures [14], inelastic scattering is not [15], and the calculated mobility in one-dimensional 
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structures in this temperature range is less than that in bulk [16]. The true origin of the 
difference lies in the fact that Dresselhaus and Rashba interactions cause a carrier's spin to 
precess slowly (during free flight) about a so-called "spin precession vector" that is defined 
by the carrier's momentum [11]. In a one-dimensional structure, a carrier is free to move 
only along one direction, and therefore the Rashba or the Dresselhaus spin precession vector 
always points along one particular direction. Scattering can change its magnitude, but not 
its direction. This leads to slow spin relaxation. In contrast, scattering can change both 
the magnitude and the direction of the spin precession vector in two- or three-dimensional 
structures. Therefore, spin depolarizes much faster in multi-dimensional structures. 

Before concluding this Letter, we should mention that in the type of structures considered 
here, there is always a magnetic field in the channel caused by the ferromagnetic contacts. 
This field, however weak, ensures that the eigenstates in the channel are not spin eigenstates 
[17]. Therefore, even non-magnetic scatterers can cause spin relaxation [18]. This mechanism 
has not been considered here, since we have not considered the channel magnetic field. In the 
absence of this mechanism, we have shown that spin relaxation length of carriers is very large 
in quasi one-dimensional structures, even at elevated temperatures and high electric fields. 
Large spin relaxation lengths have been observed before in multi-dimensional structures, but 
only at low driving electric fields and low temperatures [19]. One dimensional confinement 
can extend the range to high electric fields and elevated temperatures, which are required 
for realistic device applications. 

The work of S. P and S. B. were supported by the US National Science Foundation under 
grant ECS-0196554. 
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FIGURES 

FIG. 1. Spatial decay of the normalized spin polarized current in a GaAs quantum wire channel 
of rectangular cross-section 30 nm x 4 nm. The results are shown for four different channel electric 
fields of 1, 2, 4 and 10 kV/cm at an electron temperature of 30 K. The top panel shows schematic 
of a spin valve with a quasi one-dimensional channel. The half-metallic ferromagnetic source and 
drain contacts act as spin polarizers and analyzers, while the gate terminal is used to apply a 
transverse electric field on the channel to induce a Rashba effect. 

FIG. 2. Spatial decay of the normalized spin polarized current and the injected spin vector in 
the channel of Fig. 1. The results are shown for two different temperatures of 30 K and 77 K at a 
channel electric field of 2 kV/cm. 
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TABLES 

TABLE I. Spin relaxation length dependence on temperature and driving electric field. 



Electric field (kV/cm) 


Temperature (K) 


Spin relaxation (/xmj 


1.0 


30 


20 


2.0 


30 


60 


4.0 


30 


100 


10.0 


30 


50 


2.0 


77 


30 
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position x (|im) 




position (|im) 



